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SECTION  I 
INTRODUCTION 

The  primary  purpose  of  this  report  is  twofold.  First  it  demonstrates  the 
performance  to  be  expected  of  the  Tactical  Transmission  System  (TATS)  wideband 
modem  in  scintillated  environments  in  both  nucleaf  and  ambient  environments. 
Second,  this  study  demonstrates  some  of  the  basic  characteristics  of  UHF 
scintillated  signal  structures  as  they  are  presently  understood. 

This  report  consists  of  five  sections,  which  discuss,  in  turn,  the  environ- 
ment and  propagation,  the  TATS  system,  the  hardware  model,  the  simulation  re- 
sults, and  the  conclusions. 
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SECTION  II 

ENVIRONMENT  AND  PROPAGATION 

The  propagation  environment  geometry  assumed  here  is  representative  of  a 
link  from  equatorial  synchronous  orbit  to  midlatitude.  The  environment  is 
defined  by  the  ionospheric  mean  electron  density,  ri,  the  relative  mean  square 
electron  density  fluctuation  power,  (n-rT)^/rT^ , and  the  correlation  function  of 
the  electron  density  fluctuations.  Figure  1 contains  the  mean  electron  density 
profile  used  in  this  study  as  a function  of  both  the  altitude  and  the  line  of 
sight  (LOS)  distance  from  the  ground.  The  relative  fluctuation  power  is  varied 
to  simulate  varying  degrees  of  scintillation  effects.  The  electron  density 
fluctuation  correlation  function  is 


R(r,z) 


(n-n)^ 

g-c/L  g-(r^+a^z^)^'^^/L  _g-c/Z  g-(r^+a^z^)^/^/Z  j 

-c/L  „-c/Z 

- 

e -e 

J 

(1) 


where 


c = ZL  In  (Z/L)/(Z-L) 

Z = inner  scale  size  perpendicular  to  LOS 
L = outer  scale  size  perpendicular  to  LOS 
a = L/L' 

L‘  = outer  scale  size  parallel  to  LOS 
z = coordinate  along  the  propagation  LOS 

r = cylindrically  symmetric  coordinate  perpendicular  to  LOS 

(n-n)^  = mean  square  electron  density  fluctuation 

This  function  has  the  desirable  properties  of  having  an  inner  scale  size,  an 
outer  scale  size,  and  a zero  slope  at  r*z“0.  This  correlation  function  corre 
sponds  to  a K"^  spectrum  for  sizes  smaller  than  L but  larger  than  Z and  K"** 
for  sizes  smaller  than  z.  This  spectrum  is  characteristic  of  most  electron 
density  measurements  to  date  (refs.  1 and  2).  Sagalyn  (ref.  1)  reports  that 
this  spectrum  is  typical  of  over  90  percent  of  the  spectra  analyzed  at  AFGL. 
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This  spectrum  has  also  been  calculated  at  NRL*  in  one  of  the  most  detailed  cal- 
culations to  date. 

The  outer  scale  size,  L,  is  set  at  1.25  km,  a value  measured  by  Kelley  and 

Mozer**  in  the  equatorial  F region.  The  ratio  of  the  perpendicular  to  the 

parallel  scale  size  relative  to  the  LOS  is  assumed  as  0.36  which  is  derived  by 
considering  the  angle  (-21°)  between  the  LOS  and  the  magnetic  field  lines  in 
the  chosen  scenario.  The  inner  scale  size  is  unknown.  All  that  can  be  said 

is  that  measurements  down  to  about  30  meters  have  shown  no  evidence  of  an  inner 

scale.  For  this  calculation,  i was  arbitrarily  assumed  to  be  10  meters.  For- 
tunately, I,  as  long  as  it  is  more  than  a factor  of  10  or  more  smaller  than  L, 
has  no  significant  impact  on  the  propagation  and  subsequent  systems  studies 
assuming  that  the  spectrum  always  drops  off  K"^  or  faster.  The  assumption  that 
the  spectrum  is  monotoni cal ly  decreasing  has  been  questioned,  but  there  is  no 
evidence  that  the  higher  frequencies  contribute  significantly  in  severely 
scintillated  environments.  From  these  environmental  parameters,  the  statistics 
of  signals  propagated  through  the  environment  can  be  calculated.  Let  the  spa- 
tial dependence  of  the  signal  be 

E(z,r)  = [f(z)  + Ep(z,r)  + i Ej(z,r)j  e^*^^  (2) 


where 

z = coordinate  along  propagation  LOS 
F = vector  coordinate  transverse  to  LOS 
E = mean  field 

E[^  = in  phase  quadrature  component 
Ej  = out  of  phase  quadrature  component 

The  quantities  of  most  interest  are  the  correlation  function  of  the  quadrature 
components  which  can  be  written  as 


Eq(z,rj)  = H,^(z,r) 


(3a) 


* Scannapieco,  A.J.,  Ossakow,  S.L.,  Goldman,  S.R.,  and  Pierre,  J.M. , 
"Late-Time  Striation  Spectra,"  Naval  Research  Laboratory,  Private 
Communication,  February  1976. 

**Kelley  and  Mozer,  Private  Communication. 
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Ej(z,rJ 

Ej(z,rJ  = Hj(z,r) 

(3b) 

Er(z,FJ 

(3c) 

where  r = |r^  * >^2 1 since  the  correlation  function  of  the  electron  density 
fluctuations  is  cyl indrical ly  symmetric. 

The  F and  G functions  are  then  defined 

G(z,r)  = eTzP  + Hj^(z.r)  + Hj(z,r)  (4a) 

H(z,r)  = E(z)^  + H^(z,r)  - Hj(z,r)  + 2i  Hpj(z.r)  (4b) 

The  calculation  of  the  F and  G functions  from  the  electron  density  correlation 
function  has  been  described  elsewhere  (ref.  3)  and  will  not  be  discussed  here. 

The  decorrelation  distance  (r^)  of  the  signal  structure  is  defined  to  be  the 
1/e  point  of  G(z,r)  - E(z)^  where  z is  the  receiver  coordinate.  An  important 
limit  occurs  when  E(z)  gets  small.  Then 

Hp(z,r)  = Hj(z,r)  (5a) 

H^j(z,r)  = 0 (5b) 

Since  we  will  assume  that  Ep(z,F)  and  Ej(z,F)  are  zero  mean  gaussian  variables, 
the  above  limit  is  equivalent  to  Rayleigh  statistics.  Given  G(z,r),  F(z,r), 

E(z) , and  the  gaussian  assumption  on  the  quadrature  components,  we  have  a com- 
plete statistical  description  of  the  spatial  signal  structure. 

To  convert  spatial  statistics  to  time  statistics,  assume  that  the  character- 
istic velocity  of  the  LOS  relative  to  the  electron  density  fluctuations  is 
7 = 400  meters/second.  This  velocity  is  a reasonable  upper  limit  for  ambient 
environments  but  may  be  somewhat  low  for  certain  nuclear  environments.  The  ! 

space-velocity  transformation  is 

_ _r_ 

^7  (6)  ; 
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Procedures  for  generating  signal  samples  from  the  time  statistics  are  described 
in  the  appendix. 

Five  propagation  calculations  were  run  varying  the  relative  fluctuation 
power,  (n-ir)‘^/lT^ . Table  1 summarizes  the  calculated  environments  along  with 
three  extrapolated  environments. 

Table  1 


PROPAGATION 

CALCULATIONS 

r 

(n-n)^/rP 

E(z)^ 

To(sec) 

/ (n-n)"d 

10-" 

0.74 

2.6 

4.21x10“ 

10-' 

0.05 

1.9 

4.21x10“ 

10'^ 

4.1xl0-‘- 

0.4 

4.21x10“ 

10"' 

1.34x10-“" 

0.1 

4.21x10“ 

10° 

-0 

2.6x10"^ 

4.21x10“ 

3.07 

~0 

1.3x10'" 

1.29x10“ 

9.45 

~0 

6.5x10'' 

3.98x10“ 

29.04 

-0 

3.2x10'' 

1.22x10"’ 

where  = r^/v  is  the  decorrelation  time  of  the  signal  structure.  The  last 
three  entries  are  cases  extrapolated  from  the  = 0.026  case.  To  perform 

the  extrapolation,  (n-n)^  dz  corresponding  to  the  new  assumed  values  of 
must  be  obtained.  The  extrapolated  values  can  be  found  for  / (n-F)^dz  by 
looking  at  figure  2 which  relates  the  pertinent  quantities.  The  results  in 
this  figure  are  from  a large  number  of  calculations  varying  the  various  environ 
mental  parameters.  For  scattering  layers  where  the  environmental  parameters 
are  constant,  figure  2 can  be  simply  suimiarized. 

2.4  X 10^“  L ■ 

1 21  z L') 

V . (6) 

where  z is  the  layer  thickness  and  K is  the  signal  wave  number. 

This  extrapolation  to  different  environments  is  valid  since  the  correlation 
functions  of  the  quadrature  signal  components  of  the  = 0.026  second  case  are 
near  gaussian  in  shape.  The  gaussian  shape  is  the  limiting  form  as  the 
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propagation  effects  become  severe.  The  limiting  process  can  be  seen  by 
examining  figures  3 and  4 which  contain  the  quadrature  correlation  functions 
and  the  respective  power  spectral  densities  for  the  = 0.026  and  0.4  second 
cases.  The  evolution  toward  a gaussian  shape  is  clearly  evident.  Figure  5 
shows  the  respective  amplitude  power  spectral  densities  for  comparison.  Later, 
it  will  be  shown  by  means  of  the  system  performance  that  the  x^  = 0.4  case  is 
also  close  enough  to  the  limiting  form  to  use  for  extrapolation.  In  fact, 
all  of  the  system  results  from  all  of  the  environments  except  the  first  could 
have  been  obtained  by  extrapolation  from  any  one  of  the  environments  where 
Xq  ^ 2.  This  extrapolation,  of  course,  depends  on  the  assumption  of  gaussian 
statistics  for  the  quadrature  components.  This  assumption  is  supported  by  other 
studies  (refs.  4 and  5).  In  addition,  amplitude  spectra  of  figure  5 agree 
qualitatively  with  recent  spectral  measurements  by  Air  Force  Geophysics 
Laboratory  (AFGL)*,  adding  additional  support  to  the  Rayleigh  model. 

The  first  environment  (x^  = 2.6)  is  characteristic  of  the  rapid  transition 
from  no  effects  to  Rayleigh  statistics.  The  statistics  of  this  region  are 
not  simply  characterized  as  are  the  Rayleigh  environments.  It  does  not  matter, 
however,  as  any  successful  system  must  be  designed  to  operate  in  the  more 
stressful  Rayleigh  environments. 

The  environment  corresponding  to  x^  = 6.5  x 10"^  can  be  called  a severe 
worst  case  ambient  environment.  It  represents  a scattering  layer  400  km  thick 
with  n"  = 10®  electrons/cm^  and  (n-n’)^/^?  = 1.0.  This  case  is  also  more  severe 
than  certain  classes  of  nuclear  environments.  Thus,  the  calculations  and  re- 
sults presented  here  represent  a meaningful  range  of  real  environments. 

In  summary,  it  appears  that  severely  disturbed  UHF  signal  statistics  can  be 
simply  characterized  by  specifying  x^,  assuming  Rayleigh  statistics,  and  by 
assuming  a gaussian  correlation  function  for  the  signal  quadrature  components 
with  Xq  defining  the  1/e  point,  x^,  in  turn,  is  simply  related  to  the  basic 
environmental  parameters  via  figure  2.  This  figure,  of  course,  assumes  a K"^ 
power  spectrum.  This  simple,  but  valid,  picture  of  the  statistical  effects  of 
striated  environments  on  UHF  signals  may  perhaps  be  the  most  important  result 
in  this  report. 


*Whitney,  H. , Private  communication. 
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Figure  3.  Quadrature  Component  Correlation  Functions 
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SECTION  III 

TATS  WIDEBAND  SYSTEM  DESCRIPTION 

The  TATS  wideband  modem  is  based  on  an  8-ary  frequency  shift  key  link.  It 
is  assumed  that  the  reader  is  familiar  with  this  type  of  system.  In  addition, 
the  system  uses  a Reed-Solomon  code  to  encode  data  and  to  select  the  hopping 
of  the  carrier  frequency.  These  features  allow  multiple  user  access  via  code 
division  multiplexing  (CDM)  and  an  increased  resistance  to  noise  and  interfer- 
ence effects  such  as  jamming.  The  remainder  of  this  section  discusses  the 
transmitted  signal  structure,  the  receiver  acquisition  process,  and  the  decoding 
process  for  the  TATS  modem. 

The  signal  structure  has  two  main  features;  the  carrier  hopping  and  the 
coding-modulation  of  the  carrier.  Each  of  these  features  is  treated  in  turn 
beginning  with  the  carrier  hopping.  The  following  treatment  considers  only 
those  details  pertinent  to  the  scintillation  effects.  See  reference  6 for  a 
more  detailed  system  description. 

The  carrier  hopping  is  over  a 500  kHz  bandwidth  centered  at  260  MHz  with  the 
carrier  frequency  changing  every  10  milliseconds.  The  frequency  pattern  is 
repeated  every  seven  chips  or  every  70  milliseconds.  Each  user  has  his  own 
pattern,  the  patterns  being  chosen  to  utilize  the  500  kHz  bandwidth  and  to 
minimize  the  use  of  the  same  frequency  during  the  same  chip  by  different  users; 
that  is,  the  patterns  are  chosen  to  be  as  nearly  orthogonal  as  possible.  The 
amount  of  interference  between  users  is  determined  by  the  number  of  users  using 
the  system  bandwidth  at  any  one  time.  This  use  of  CDM  has  the  advantage  that 
as  the  number  of  users  decreases,  the  link  margin  of  the  remaining  users  in- 
creases insuring  efficient  use  of  the  system  power  and  bandwidth. 

It  has  been  shown  elsewhere  (ref.  7 and  NELC*  data)  that  coher- 
ence bandwidths  in  frequency  space  for  scintillations  are  much  wider  than 
500  kHz  so  it  is  reasonable  to  ignore  frequency  selectivity  over  the  hopping 
bandwidth.  The  degradation  of  the  signals  due  to  the  presence  of  other  users 
is  also  ignored  in  this  study.  These  two  assumptions  mean  that  the  details  of 

♦Private  communication. 
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the  carrier  frequency  hopping  need  not  be  explicitly  modeled.  The  details  of 
the  information  coding,  however,  must  be  modeled  in  detail. 

The  information  coming  into  the  modulator  is  in  terms  of  6 bit  characters. 
These  characters  are  encoded  using  a Reed-Sol omon  code  where  each  of  the  6 bit 
characters  is  represented  by  a sequence  of  seven  numbers,  each  number  ranging 
from  0 to  7.  The  sequences  for  each  character  are  chosen  to  minimize  the 
occurrence  of  the  same  number  in  the  same  position  in  the  7 number  sequence 
for  any  two  characters,  that  is,  maximize  character  orthogonality.  For  the 
code  used  here,  a particular  character  representation  has  no  numbers  in  its 
sequence  in  common  with  14  other  characters  and  only  one  number  in  common  with 
the  remaining  49  characters.  The  use  of  these  nearly  orthogonal  representations 
results  in  a significant  improvement  in  performance  over  noncoded  schemes. 

A coded  character  is  modulated  onto  the  carrier  frequency  in  a sequence 
of  eight  chips,  the  last  seven  chips  corresponding  to  the  sequence  of  seven 
numbers  of  the  Reed-Solomon  character  representation.  The  first  chip  is  at 
the  carrier  frequency  and  is  used  to  track  time  and  frequency  changes.  Each 
chip  is  modulated  using  an  8-ary  frequency  shift  key  (FSK)  modulation.  In 
this  modulation,  each  of  the  eight  possible  number  values  in  the  Reed-Solomon 
code  is  represented  by  a discrete  frequency  relative  to  the  carrier  frequency. 
For  example,  the  6 bit  character  100101  is  represented  by  the  sequence 
(1 ,2, 3, 7, 5, 0,6).  The  code  word  representation  in  frequency  space  is  shown 
in  figure  6.  "f^,"  is  the  carrier  frequency  and  Af  = 100  Hz. 

It  is  important  to  note  that  a code  word  is  one  chip  longer  than  the 
carrier  frequency  hopping  pattern.  As  a result,  the  synchronization  chip  at 
the  carrier  frequency  is  rotated  through  the  carrier  hopping  pattern  giving 
a measure  of  frequency  diversity  protection  to  the  synchronization  process. 

The  signal  structure  development  for  information  transmission  can  be 
summarized  as  follows.  During  transmission,  a 6 bit  information  character  is 
put  into  the  coder,  which  converts  it,  using  a Reed-Solomon  code,  to  a sequence 
of  seven  numbers  ranging  from  0 to  7.  The  Reed-Solomon  sequence  is  modulated 
on  seven  chips  of  the  carrier  frequency  using  an  8-ary  FSK  modulation.  An 
unmodulated  synchronization  chip  precedes  each  7 chip  sequence  and  the  8 chips 
together  represent  the  original  6 bit  information  character.  Finally,  the 
carrier  frequency  is  hopped  over  a 500  kHz  bandwidth  every  chip.  During 
acquisition,  the  signal  structure  is  somewhat  different.  First  there  are  532 
chips  of  unmodulated  hopped  carrier.  Since  each  user  has  his  own  pattern,  this 
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chip  sequence  acts  as  a user  address.  The  user  employs  the  sequence  to  acquire 
and  pull  in  both  time  and  frequency  synchronization.  Four  start-of-message 
characters  follow  this  chip  sequence.  These  characters  have  the  same  structure 
as  information  characters  except  no  synchronization  chips  are  used.  These 
characters  are  also  sent  having  each  7 chip  sequence  in  time  with  the  7 chip 
carrier  frequency  hopping  pattern. 

Next,  it  is  necessary  to  know  the  procedure  used  by  the  demodulator  to 
recognize  the  signal.  The  hopping  pattern  recognition  process  is  carried  out 
by  eight  matched  filters  that  are  also  used  to  decode  the  information  charac- 
ters. The  filters  are  separated  by  100  Hz;  thus,  the  frequency  uncertainty 
that  can  be  initially  handled  is  about  ±400  Hz.  (In  practice,  this  range  can 
be  increased  by  shifting  the  initial  carrier  frequency  estimate  to  widen  the 
coverage  of  the  eight  matched  filters  in  frequency  space.)  When  searching 
for  its  pattern,  the  receiver  dehops  the  carrier  and  feeds  the  result  to  eight 
matched  filters  each  of  which  is  followed  by  an  envelope  detector.  The  quan- 
tized outputs  of  each  filter-envelope  detector  are  sampled  at  T/2  and  T relative 
to  the  receiver.  The  outputs  are  accumulated  in  16  accumulators,  2 for  each 
filter.  If,  at  the  end  of  21  chips,  none  of  the  16  accumulators  has  exceeded 
a preset  threshold  of  110,  the  accumulators  are  cleared,  the  dehopping  pattern 
is  slipped  a chip  in  timing,  and  another  21  chip  search  is  done.  This  process 
continues  until  the  local  dehopping  pattern  is  within  about  T/2  seconds  time 
synchronization  with  the  sender  hopping  pattern.  When  this  occurs,  given 
sufficient  power  and  given  that  the  carrier  frequency  is  not  in  error  by  more 
than  about  ±400  Hz,  one  of  the  16  accumulators  should  exceed  the  preset  thresh- 
old because  at  least  one  of  the  filters  must  be  close  to  the  frequency  being 
sent  and  the  filter  integration  ending  at  T/2  or  T must  overlap  the  signal  sent. 
When  an  accumulator  exceeds  the  threshold,  the  accumulator  address  is  stored, 
the  accumulators  are  cleared,  and  another  21  chip  search  is  initiated  but 
without  changing  the  dehopping  timing.  If  the  same  accumulator  again  exceeds 
the  threshold  within  21  chips,  the  receiver  assumes  its  pattern  is  being  sent. 
The  532  chips  allow  for  two  chances  to  recognize  the  pattern  before  a failure 
will  occur.  If  the  pattern  recognition  occurs  near  the  end  of  the  532  chips, 
sufficient  chips  may  not  remain  to  finish  the  acquisition  process. 

Once  the  pattern  has  been  recognized,  the  receiver  does  an  initial  time 
and  frequency  correction  based  on  the  address  of  the  accumulator  that  exceeded 
the  threshold.  If  the  accumulator  corresponded  to  a sample  at  T,  the  timing 
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is  retarded  by  T/4;  otherwise  no  correction  is  made.  The  frequency  shift  of 
the  filter  where  the  signal  was  found  is  added  to  the  initial  carrier  frequency 
estimate  correcting  that  quantity.  With  these  initial  corrections,  the  re- 
ceiver proceeds  to  pull  in  the  time  and  frequency  track  more  accurately. 

The  time  and  frequency  tracking  are  done  with  three  matched  filters.  For 
frequency  tracking  two  filters,  tuned  to  f^  tiif/2  where  f^  is  the  carrier 
frequency  estimate  resulting  from  the  hopping  pattern  recognition  process, 
are  used.  These  filters  integrate  over  the  chip  period  and  the  quantized  out- 
puts are  compared.  If  the  f + Af/2  output  is  greater,  a counter  is  increased 

C i 

by  one.  If  the  f^  - Af/2  output  is  greater,  the  counter  is  decreased  by  one. 

If  the  counter  value  becomes  greater  than  NC,  the  counter  is  zeroed  and  f^  is 
increased  by  Af/8.  If  the  counter  value  is  less  than  -NC,  f^  is  corrected  the 
other  way.  Thus,  the  carrier  frequency  estimate  is  continually  adjusted 
toward  the  real  carrier  frequency.  The  timing  is  handled  much  the  same  way. 

A filter  tuned  to  the  carrier  frequency  estimate  is  sampled  at  T/2  and  T and 
the  quantized  outputs  are  compared.  If  the  T output  is  greater,  another 
counter  is  increased  by  one.  If  the  T/2  output  is  greater,  the  counter  is 
decreased.  The  timing  is  corrected  by  ±T/32  whenever  the  timing  counter 
exceeds  NC  or  is  less  than  -NC,  respectively.  In  this  tracking  mode,  the 
system  tracks  for  98  chips  with  NC  = 3 and  then  for  another  98  chips  with 
NC  = 7.  After  the  196  chips  have  been  processed,  the  timing  and  frequency 
are  presumably  sufficiently  accurate  to  look  for  the  start  of  message  charac- 
ters. 

The  system  looks  for  the  start  of  message  characters  for  343  chips  and 
goes  back  to  the  hopping  pattern  search  if  they  are  not  found.  The  search  for 

the  4 characters  is  done  by  the  same  procedure  as  the  decoding  of  the  informa-  : 

tion  characters  with  two  differences.  First,  there  are  no  synchronization  I 

chips.  Second,  the  criterion  for  acceptance  includes  the  sum  of  the  4 start-of-  ] 

message  character  accumulators  being  greater  than  245  as  well  as  recognition  ' 

of  the  individual  characters.  After  these  characters  are  found,  the  system  ’ 

has  time,  frequency,  and  character  synchronization,  and  information  decoding 

can  begin.  i 

When  information  is  being  received,  the  receiver  is  carrying  out  three 
basic  functions.  First,  it  maintains  time  and  frequency  synchronization  using 
the  synchronization  chip  at  the  beginning  of  each  character.  The  procedure  for 
this  chip  is  identical  to  that  used  in  the  initial  pull-in  during  acquisition 
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with  NC  = 7.  The  maximum  tracking  rates  are  easy  to  estimate.  The  timing  loop 
can  track  up  to  about  5xl0‘‘*  sec/sec.  The  frequency  loop  can  track  up  to  about 
19.5  Hz/sec. 

The  second  and  most  important  function  is  to  decode  the  Reed-Solomon  code. 
To  each  of  the  64  characters,  there  corresponds  an  accumulator  represented 
schematically  in  figure  7.  As  each  chip  is  processed  through  the  matched 
filter  bank,  each  character  accumulator  receives  the  quantized  output  from  the 
appropriate  filter.  For  example,  the  100101  character  discussed  earlier  is 
the  38th  character;  thus,  the  38th  accumulator  receives  successively  quantized 
outputs  from  filters  1,2, 3, 7, 5,0,  and  6.  At  the  end  of  the  seven  chips,  all 
of  the  accumulators  are  compared  and  the  character  corresponding  to  the  accumu- 
lator with  the  largest  value  is  inputted  into  the  data  stream.  For  the  present 
work,  random  characters  were  used  and  the  results  are  in  terms  of  6 bit  charac- 
ter error  rates. 

The  final  function  performed  during  information  decoding  is  detecting 
loss  of  message.  This  is  accomplished  by  sampling  the  largest  value  from 
the  character  accumulators  at  the  end  of  each  character.  If  three  out  of 
four  successive  values,  where  one  of  the  three  must  be  the  first,  is  below 
a threshold  of  70,  the  system  starts  a time-out  counter.  If  four  successive 
character  values  that  exceed  the  threshold  are  not  found  within  6 seconds, 
the  system  stops  decoding  and  returns  to  a pattern  hopping  search  mode.  If 
the  four  character  values  are  found,  the  system  returns  to  its  normal  decoding 
state. 
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Figure  7.  8-ary  Demodulator 
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SECTION  IV 
HARDWARE  MODEL 

As  shown  in  the  previous  section,  the  acquisition,  decoding,  and  other 
system  functions  are  based  on  comparisons  of  the  quantized  outputs  of  the 
matched  filter-envelope  detector  hardware.  To  calculate  these  outputs  it  was 
necessary  to  develop  computer  algorithms  simulating  the  system  hardware 
including  the  automatic  gain  control  (AGC).  This  section  considers  the  model 
formalism  used  to  provide  an  accurate,  efficient  simulation  of  the  hardware 
functions.  First,  the  matched  filter  model  is  described  as  follows. 

The  received  signal  is  represented  as 

S(t)  = E(t)e^'^^^^c  + n(|.)  (7 

where 

E(t)  = I + Ep(t)  + i Ej(t) 
f^  = carrier  frequency 
i = il , ±3,  i5,  ±7  = frequency  shift  index 
/if  = filter  frequency  separation 
n(t)  = white  gaussian  noise  input 

E(t)  is  the  same  as  E(r,z)e"  in  equation  2 with  t = |r|/v  and  z equal  to 
the  receiver  coordinate.  The  output  of  the  matched  filter  corresponding  to 
the  mth  frequency  shift  with  the  nth  shift  present  is 


R 


mt 


+ fg]t 


+ 


. f 


(8) 


21 


AFWL-TR-76-128 


where 

N = unilateral  power  spectral  noise  density 
= undisturbed  field  level 
= demodulator  carrier  frequency  estimate 

% - 'c-'c 

T'  =T-iTel 

I = chip  period 
Tg  = timing  error 

The  filter  output  in  equation  9 is  divided  into  three  contributions.  The 
first  two  result  from  dividing  the  first  integral  into  low  and  high  frequency 
signal  portions  and  the  last  is  noise.  The  low  frequency  portion  of  E(t)  is 
generated  using  Fourier  techniques  described  in  the  appendix  and  is  represented 
as  a Fourier  series. 


,nTTt 


n=-n. 


(9) 


where  n^  = L/T ; 2L  is  the  time  domain  of  the  Fourier  expansion  and  is  chosen 
much  larger  than  the  correlation  time,  of  the  signal  structure.  The  resolu- 
tion of  Ej^p(t),  it,  is  chosen  to  be  greater  than  or  equal  to  T.  During  each 
it,  assumed  linear.  Let 

(Elf(^'*’^^)  “ E|^p(t)) 

^ " At  (10a) 


b = E^pCt) 


OOb) 


The  contribution  to  from  ^LF  (t)  can  be  written  as 
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where  the  range  of  the  integral  is  in  the  At  interval  fit  by  a and  b.  After  a 
change  of  variables. 


J‘ 


m£ 


. ^ „i2TT[(£-m)Af/2  + f. 
■ N ^ e 


LF 
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(12) 


Finally, 
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where 


G (!•  f 1 = ^^0  . /-i2Tr[(£-m)Af/2  + f ]T'  , \ 

'^e^  Mf(m-£)N  ® ^ ) 


(14) 
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irAf  ( 
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^ -l) 


(15) 


b'  = b + a t' 


(16) 


Thus,  the  low  frequency  part  of  the  filter  output  is  reduced  to  a sample  expres- 
sion. The  H and  G functions  need  be  evaluated  only  every  At  which  may  be  much 
larger  than  T.  It  will  become  clear  later  that  the  phase  of  is  not 

important.  Thus  the  phase  factor  in  equation  (14)  need  not  be  evaluated.  If 
the  filters  are  identical,  the  G and  H functions  become  only  functions  of  m-£ , 
which  then  requires  only  15  separate  H and  G coefficients  rather  than  64. 
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The  high  frequency  contribution  to  the  filter  output  is 


mil 


HF 


5C  i:  * V'dt  (17) 

■^0  n=n_+l 


These  frequencies,  since  they  are  faster  than  a chip  period  are  assumed  to 
contribute  as  random  gaussian  noise  in  the  filters.  After  integrating  and 
averaging,  the  variance  of  the  filter  output  is 
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(18) 


(TSf,)  (19) 

C »C 

where  i s the  signal  power  spectral  density  and  Sinc(x)  = sin^(x)/x^.  This 

treatment  of  the  higher  frequency  components  assumes  that  R^^ R^^  equals  zero 
and  that  the  filter  outputs  from  these  components  for  each  successive  chip  are 
independent.  The  first  assumption  is  always  true  since  propagation  studies 
demonstrate  that  the  signal  structure  is  Rayleigh  before  significant  power  is 
found  at  these  higher  frequencies.  The  second  assumption  is  an  approximate 
result  of  specifying  n^  = L/T. 

The  final  filter  contribution  is  due  to  noise,  which  is  assumed  to  be 
gaussian,  white,  and  equal  in  all  filters.  It  is  easy  to  demonstrate  that 
the  noise  output  has  a variance  of 


2d^  = 2Eq^T/N  (20) 

which  is  four  times  the  chip  energy  to  noise  ratio.  The  term  can,  for  conve- 
nience, be  added  to  Rj^*R^j^  since  both  are  variances  of  independent  gaussian 
variables. 

The  output  of  the  mth  filter  where  the  signal  was  transmitted  in  the  tth 
is 
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1 


* (b  ♦ af)  G„,(T’,f^)  ♦ (F„j(T',f^)  ♦ d=)>/MG(o)  * 1 G(o)) 

(21) 


where  G(o)  is  a gaussian  sample  of  unit  variance.  This  particular  model  allows 
for  the  effects  of  timing  errors,  frequency  errors,  and  scintillation  induced 
phase  and  amplitude  effects  without  being  complex  or  expensive  to  execute. 

The  matched  filter  outputs  are  each  envelope  detected,  which  corresponds 
to  taking  the  absolute  value  of  The  resulting  scalar  quantities  are 

quantized  to  integer  quantities  with  a maximum  of  15  and  then  processed  for 
acquisition  or  information  decoding. 

The  next  piece  of  hardware  to  be  considered  is  the  automatic  gain  control 
(AGC).  The  AGC  maintains  the  signal  power  into  the  demodulator  filters  so  that 
the  average  envelope  detection  output  after  quantization  is  5.5  with  only 
noise  present.  As  S/N  becomes  very  large,  assuming  fg  = T^  = 0,  the  output 
of  the  envelope  detector  following  the  filter  containing  the  signal  approaches 
25.  Note  that  the  quantization  truncates  this  value  to  15.  It  is  necessary 
to  model  the  AGC  because  some  of  the  information  processing  is  dependent  on 
thresholds  for  the  quantized  outputs.  The  dynamic  response  of  the  AGC  is 
mcxJeled  as  exponential.  For  a step  function  increase  in  signal  strength,  the 
gain  signal  decreases  exponentially  to  its  new  level  with  a time  constant  of 
16  milliseconds.  For  a step  decrease,  the  time  constant  is  80  milliseconds. 
This  asymmetry  in  response  degrades  the  system  when  the  signal  fluctuations 
become  comparable  to  these  response  times  since  the  signal  level  is  depressed 
by  the  asyimietry.  This  question  was  not  examined  in  detail  since,  as  will  be 
shown  later,  the  severest  system  degradation  occurs  for  circumstances  where 
AGC  effects  are  not  significant. 

The  modeling  of  the  system  processes  is  straightforward  given  the  hardware 
models  already  described.  First,  using  Fourier  techniques  found  in  the  appen- 
dix, sample  signal  structures  are  generated.  This  signal  is  then  used  to  cal- 
culate the  relevant  AGC — filter— envelope— detector—  quantizer  outputs  for 
each  successive  chip.  These  outputs  are  processed  according  to  the  description 
of  the  acquisition  and  decoding  procedures  discussed  in  the  last  section.  The 
results  are  finally  accumulated  in  terms  of  various  statistical  measures  shown 
in  section  V. 
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SECTION  V 

SIMULATION  RESULTS 

Figure  8 contains  the  results  of  the  acquisition  study.  The  results  are 
characterized  by  four  regions.  The  first  (0.7  r^/  v _<  £ 1.4  r^/  v)  is  a 

transition  region  from  no  degradation  effects  to  relatively  constant  effects. 
The  transition  region  is  relatively  small  and  is  located  at  7.  The 

signal  statistics  in  this  region  are  generally  nonRician  and  not  amenable  to 
simple  descriptions  as  are  the  other  regions. 

The  second  region  (0.2  1 _<  0.71^^/  7)  of  relatively  constant  system 

performance  is  the  most  important  region  as  it  represents  the  worst  system 
performance  over  most  of  the  likely  range  of  x^^.  The  relatively  constant 
performance  is  a result  of  the  fact  that  the  signal  fluctuations  are  slow 
compared  to  the  effective  integration  periods  of  the  acquisition  algorithms. 

The  third  region  (0.03  1 Tq  1 0*2)  shows  a recovery  of  the  acquisition 
performance  caused  by  the  integration  properties  of  the  acquisition  algorithms 
as  Xq  becomes  smaller.  Finally,  the  fourth  region  (x^  _<  0.02)  shows  a system 
degradation  caused  by  an  increase  in  scintillation  induced  cross-talk.  The 
cross-talk  becomes  important  when  l/x^  approaches  the  frequency  separation  of 
the  matched  filters. 

This  study  indicates  that  about  18  dB  of  S/N  are  required  to  get  a 0.84  or 
greater  probability  of  acquisition  for  any  ambient  environment  at  18  dB  or 
above.  Most  of  the  failures  in  acquisition  were  due  to  failure  in  recognizing 
the  start  of  message  characters.  The  initial  hopping  pattern  recognition  and 
the  time  and  frequency  pull-in,  in  general,  worked  well  even  with  a severely 
scintillated  signal. 

Figure  9 shows  the  ambient  acquisition  performance  along  with  the  perfor- 
mance in  the  x^  = 1.9  second  environment.  The  second  curve  represents  a 
limiting  case  of  the  acquisition  performance  for  any  nonnuclear  environment. 

After  successful  acquisition,  the  system  must  be  able  to  decode  the  infor- 
mation in  the  signal.  Three  aspects  of  this  problem  were  examined.  First, 
the  ability  to  maintain  time  and  frequency  tracking  was  considered.  None  of 
the  environments  considered  here  seriously  affected  the  time  and  frequency 
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tracking  behavior.  Second,  the  decoder  must  correctly  recognize  the  characters 

sent.  Figure  10  shows  the  6 bit  character  error  rate  as  a function  of  t and 

0 

S/N.  The  performance  curves  are  similar  in  character  to  the  acquisition  curves 
for  the  same  basic  reasons.  Figure  11  has  the  ambient  performance  along  with 
the  results  of  the  = 0.013  and  1.93  second  environments.  The  = 1.93 
curve  can  be  considered  as  a limiting  performance  curve  for  any  ambient  environ- 
ment. In  using  the  curves  in  figures  10  and  11,  it  must  be  remembered  that  the 
character  errors  are  not  random  in  time.  Figures  12  and  13  show  the  simulation 
results  for  a string  of  characters  at  S/N  of  18  dB  and  13  dB  with  x^  = 1.9 
seconds  where  the  dashes  represent  incorrectly  received  characters.  The  non- 
random  nature  of  the  errors  is  clearly  evident.  The  effects  of  the  nonrandom- 
ness of  the  errors  depends  on  the  nature  of  the  message  itself;  but,  in  general, 
nonrandom  errors  are  not  desirable. 

Earlier,  in  section  II,  it  was  shown  that  propagation  effects  can  be 
extrapolated  from  a single  propagation  calculation  if  that  calculation  is 
reasonably  near  the  gaussian  limit  of  the  correlation  functions  of  the  quadra- 
ture components.  Figure  10  shows,  for  S/N  = 18  dB,  the  system  performances  for 
a few  cases  extrapolated  from  the  x^  = 0.4  case.  The  extrapolated  results  agree 
well  with  the  nonextrapolated  calculations.  The  extrapolated  result  at  x^  = 4 
is  what  would  have  happened  if  7 were  a bit  smaller  or  r^  a bit  larger  causing 
the  onset  of  effects  to  occur  at  larger  x^. 

The  numerical  error  in  figure  10  and  those  following  is  primarily  due  to 
jitter  from  the  Monte  Carlo  aspect  of  the  calculation.  For  x^  ^ 0.05,  the 
estimated  error  is  less  than  20  percent  relative;  and  for  x^  < 0.05,  the 
estimated  error  is  less  than  50  percent  relative.  A factor  of  about  2 was 
added  to  these  error  estimates  over  the  actual  estimated  numerical  error  to 
account  for  the  uncertainties  in  the  numerical  modeling.  Unfortunately,  avail- 
able system  performance  information  was  insufficient  to  evaluate  the  accuracy 
of  the  numerical  models  even  with  nonperturbed  signals. 

The  third  aspect  of  the  decoding  examined  was  the  loss  of  message  algorithm. 
Figure  14  shows  the  average  turn-off  time  for  S/N  = 18  dB.  For  0.2  1 Tq  £ 0.7 
r^/V,  the  system  would  disable  itself  after  only  10  seconds  on  the  average. 

This  is  clearly  not  acceptable,  but  fortunately  the  loss  of  message  threshold, 
CTR,  is  adjustable.  Figure  15  shows  the  average  turn-off  time  as  a function 
of  CTR  with  S/N  = 18  dB  and  x = 0.4.  The  problem  can  be  corrected  by  lowering 
CTR  but  at  the  cost  of  not  recognizing  the  loss  of  signal  as  readily.  Since 
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the  lower  limit  of  CTR  is  48,  it  seems  clear  that  18  dB  of  S/N  is  required  for 
reasonable  system  performance.  Finally,  figure  16  shows  the  average  turn-off 
time  for  CTR  = 70  as  a function  of  S/N.  The  S/N  must  approach  30  dB  to  get 
messages  through  that  are  more  than  a couple  of  minutes  long  in  the  nominal 
system  configuration  used  here. 
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SECTION  VI 

CONCLUSIONS  AND  RECOMMENDATIONS 

1.  The  loss  of  message  algorithm  posed  the  most  serious  problem.  For 
the  nominal  configuration  assumed  for  this  study,  an  S/N  of  about  28  dB  is 
required  to  get  message  lengths  of  2 minutes  or  longer.  This  problem  can  be 
alleviated  by  lowering  the  loss  of  message  threshold.  With  the  minimum  setting 
of  the  threshold,  messages  of  several  minutes  duration  will  be  received  at 

18  dB.  It  is  recommended  that  this  problem  be  examined  further,  preferably  by 
testing  the  real  hardware. 

2.  Neglecting  the  loss  of  message  problem,  the  system  worked  marginally 
well  at  18  dB.  The  probability  of  acquisition  was  0.84  and  the  6 bit  character 
rate  was  about  0.05.  Operation  at  this  level  requires  toleration  of  character 
error  behavior  as  shown  in  figure  12.  Figures  9 and  11  contain  worst  case 
system  performance  curves  for  any  perturbed  nonnuclear  environment. 

3.  The  ultimate  limit  in  system  performance  as  the  environment  becomes  more 
severe  is  set  by  the  chip  period  and  the  filter  separation.  If  either  « T 
or  Tq  <<  1/Af,  then  the  system  performance  is  irrecoverable  for  practical  S/N. 

4.  The  use  of  the  Reed-Solomon  code  resulted  in  a significant  improvement 
in  system  performance  due  to  the  orthogonality  of  the  character  representations 
and  the  integration  effect  of  the  code  over  the  character  length. 

5.  UHF  propagation  is  sufficiently  well  understood  that  worst  case  per- 
formance curves  can  be  developed  for  UHF  systems.  This  means  that  system  speci- 
fications and  the  means  to  test  for  compliance  are  now  practical  for  scintil- 
lation effects. 

6.  It  is  recommended  that  a signal  processor  be  developed  to  simulate 
scintillation  effects  for  system  test  purposes.  Because  of  the  complexity  of 
these  systems,  it  is  not  likely  that  sufficient  confidence  in  survivability 
can  be  gained  from  simulations  such  as  done  here. 
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APPENDIX 

GENERATION  OF  COMPLEX  CORRELATED  RANDOM  SEQUENCES 


Let  f(t)  be  a zero  mean  random  gaussian  variable  on  the  domain  (-T,T). 
f(t)  can  be  expanded  in  a Fourier  series. 


■ 7 Z “n  ■ f 


-a; 


f(t)  dt 


G(t-f)  = f*(t)f(f) 


F(t-f)  = f(t)f(f) 


Lim(TC*  C|^)  = 0,  m^^n.m^^O 


Lim(TC*  C|^)  = 4 / da  G(a)  cos  (oaij^) 


Lim(TC^  Cjj,)  = 0,  m n,  m 0 


Lim(TC^  C_|^)  = 41  da  F(a)  cos  (aU|^) 


= T(L„  + i M„) 
' m m' 
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If  T is  chosen  such  that  F(t)  = G(t)  - 0 for  t T,  the  limits  in  equations  (A5), 
(A6),  (A7),  and  (A8)  can  be  assumed.  Now  if  + iC^^. , and  using 

equations  (A5),  (A6),  (A7),  and  (A8),  then  with  m ^ 0, 


C C 
mr  mr 


• = C^i  C^i  = V2 

(A9) 

mr  ^mi  " ^ 

(AlO) 

= -^mi 

(All) 

= '-mr  'mi  = 

(A12) 

Since  and  are  random  gaussian  variables,  the  joint 

probability  distribution  function  of  these  variables  can  be  written  in  terms 

find  M^. 
m m’  m 

mi)  " ^ 6XP 

mr  mi  -mr  -mi' 


Km 

0 

-'m 

-Mm 

'mr 

ic  c . c c . 1 

rmr  mi  -mr  -mi ' 

0 

Km 

-Mm 

'mi 

-'m 

-Mm 

Km 

0 

'-mr 

-Mm 

'm 

0 

Km 

'-mi 

(A13) 


This  expression  can  be  simplified  by  diagonalizing  the  4x4  matrix.  To 
diagonalize  the  above  matrix,  a transformation  matrix  must  be  found  with  the 
following  properties. 


T-iT  = jj-i  = I 

T Q T"^  = R 
40 


(A14) 

(A15) 


I 


AFWL-TR-76-128 


where  I is  the  identity  matrix,  Q is  the  4x4  matrix  in  equation  13,  and  R is 
a diagonal  matrix  with  diagonal  elements  X^,  X^,  and  X^.  The  diagonal 
elements  are  easily  solved  for 


X , = X = 1C  + s/L^  + M 
1 2 m m I 


^3  = ^ 

A matrix,  T,  with  the  desired  properties  is 


where  sf 
as 


T = 


VT 
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(A18) 


probability  density  function  can  now  be  written 


where 
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(A19) 


(A20) 

The  primed  variables  are  also  random  gaussian  variables  with  variances  as  shown 
in  equation  (A19).  Generating  a simple  of  f(t)  on  the  interval  -T  £ t £ T is 
now  straightforward.  For  each  m 0,  sample  a gaussian  distribution  four  times 
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to  get  values  of  C V , and  , transform  using  equation  (20)  to 

unprimed  variables,  and  form  the  random  Fourier  coefficient  by  C = C + iC  .. 

m mr  mi 


For  m = 0,  another  procedure  must  be  used.  From  equations  (A6)  and  (A8J 


C"  + 
or  01 


C"  - 
or  01 


= L, 
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C C . = M /2 
or  01  o' 
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Thus 
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(A24) 
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The  joint  probability  function  can  now  be  written  as 
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(A-26) 


Following  the  same  procedure  as  with  m ^ 0,  the  probability  matrix  can  be 
diagonal ized. 


Let 


[(Kq+Lq)/2]i/^ 

K * 

(A-27) 

[(Kq-Lq)/2]‘/^ 

(A-28) 
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Then 


(A30) 


To  generate  and  , a random  gaussian  sample  of  and  are  generated 
with  the  variances  and  respectively,  and  used  in  equations  (A27)  and 
(A28).  After  the  random  complex  Fourier  coefficients  are  found  for  all  m,  the 
sample  f(t)  is  calculated  by  inverse  Fourier  transforming  as  in  equation  (Al). 
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